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From a mean-field solution of the Hubbard-Holstein model, we show that a rich variety of different 
electronic phases can result at the interface between two polar materials such as LaAlOs/SrTiOs. 
Depending on the strengths of the various competing interactions, viz., the electronic kinetic energy, 
electron-phonon interaction, Coulomb energy, and electronic screening strength, the electrons could 
(i) either be strongly confined to the interface forming a 2D metallic or an insulating phase, (ii) 
spread deeper into the bulk making a 3D phase, or (iii) become localized at individual sites forming 
a Jahn- Teller polaronic phase. In the polaronic phase, the Coulomb interaction could lead to 
unpaired electrons resulting in magnetic Kondo centers. Under appropriate conditions, electronic 
phase separation may also occur resulting in the coexistence of metallic and insulating regions at 
the interface. 

PACS numbers: 73.20.-r, 71.38.-k, 64.75.Jk 
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I. INTRODUCTION 

Transport measurements of the polar interfaces such as 
LaAlOg/SrTiOg (LAO/STO) have revealed a rich variety 
of behavior^ such as a Kondo resistance minimum, metal- 
lic, insulating, and superconducting characters, as well as 
magnetism under varying experimental conditions. Su- 
perconductivity has also been observed^ below 200 mK 
and while the resistance minimum can be explained by 
either the 2D weak localization or the presence of Kondo 
centers, it has been argued that the observation of a neg- 
ative magnetoresistance together with the logarithmic in- 
crease of resistivity with decreasing temperature suggests 
the presence of localized Kondo centersi^ The resistance 
also follows the scaling laws for the Kosterlitz-Thouless 
transition in 2D," which can be tuned by the carrier con- 
centration induced by a gate voltage.^ Yet another in- 
triguing experimental discovery is the reversible switch- 
ing between nanoscale conducting and insulating regions 
and patterning of conducting lines on the surface of the 
polar heterostructure using an atomic force microscopy 
probC)^ 

The origin of the 2D electrons at the polar interfaces is 
becoming increasingly clear to be due to electronic recon- 
struction necessary to quench the polar catastrophe, i.e., 
the diverging electrostatic potential due to the lining up 
of the charged atomic planes must be mitigated by trans- 
ferring electrons to the interface. However, the complex 
behavior of these electrons originating from the compe- 
tition between the various interactions such as Coulomb 
and Jahn- Teller (JT) terms and external factors such as 
impurities and gate voltage that control the carrier den- 
sity is less understood. 

In this paper, we study a model polar system using 
mean-field theory and examine the various phases of the 
electrons at the interface. We find that depending on 
the strength of the various competing interactions, the 
interfacial electrons show a rich variety of phases. Under 
appropriate circumstances, a phase separation between 



metallic and insulating phases is also possible. 



II. MODEL AND METHOD 

Before discussing the model Hamiltonian, we present 
a brief overview of the relevant electron states at the po- 
lar interfaces, taking the example of LAO/STO for the 
sake of concreteness. Simple electron counting argument 
shows that half an electron per interfacial Ti atom are 
needed for the ideal interface to quench the divergence of 
the electrostatic potential. The electrons form a quasi- 
2D electron gas on the STO side, where they occupy 
the Ti(d) bands as indicated from density-functional 
calculations^"— as well as a variety of experimental. 
They furthermore become coupled to the lattice, as the 
Ti in SrTiOs, originally in the d° state, becomes a d^ 
JT ion. Crystal field splits the Ti t2g states below 
which then interact via the JT interaction with several 
octahedral stretching and shear modes, leading to a com- 
plex distortion of the octahedra at the interface. In the 
present work, we restrict to the simplified Holstein model 
of the electron-lattice interaction, which describes a sin- 
gle orbital degree of freedom. The model captures the 
essential physics of the competing interactions including 
the polaron formation, although effects such as orbital 
ordering are obviously not contained in the model. 

The Hubbard-Holstein model adopted here describes 
the motion of the electrons on the Ti sites with the other 
ionic charges serving as a Coulombic background. The 
Hamiltonian is 

Ti. = TiMad + ^^e + He-p, (1) 

where the first term describes the Madelung energy 



ij Qi Qj , 



(2) 



M being the Madelung matrix, which may be calculated 
using the standard Ewald's method, i and j are the site 
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FIG. 1: Model of a typical polar interface considered in the 
paper (Top). Here we are concerned with the n-type structure 
with Ti02 and LaO forming the interfacial layers. Electrons 
migrate to the interface in response to the diverging Coulomb 
potential of the charged planes and interact via Coulomb 
and JT interactions, resulting in a rich variety of interfacial 
phases. Bottom part shows the calculated bare and screened 
Coulomb potentials and the accumulated electron density at 
the interface in response to the polar catastrophe (indicated 
by the bare potential, dashed line). Potentials shown are the 
computed Madelung potentials on the oxygen sites in the SrO 
or LaO layers. 



indices, qi is the total of the ionic charge and the elec- 
tronic charge rii on the i-th site (qi = q'^ + rii), e is the 
static dielectric constant, and the prime on the summa- 
tion excludes same site terms. Note that while the site 
summations run over all atoms, the electronic charge rii 
is non-zero on the Ti sites only. Without the electronic 
charges the Madelung energy would diverge as the po- 
tential grows unrestricted away from the interface (polar 
catastrophe). The various ions are passive, merely pro- 
viding a Coulomb potential background for the Ti elec- 
trons to move. 

The second term represents the electronic kinetic en- 
ergy and the electron-electron interaction 

He= c\^Cj„ + H.c. + ^ ni^rii^, (3) 

where Cia-, cj^ are the field operators, a is the spin in- 
dex, and {ij) indicates summation over distinct nearest- 
neighbor pairs. An on-site Coulomb repulsion (U) is re- 
tained in the electron-electron interaction term. The final 
term in the Hamiltonian is the electron-lattice interaction 



within the Holstein model 

'ne-p^Y.^\KQj~gQ,n,), (4) 

i 

where Q is the lattice distortion at a lattice site, K and 
g are the stiffness constant and the electron-lattice cou- 
pling strength, respectively, and rii is the electron occu- 
pancy of the Ti atom. In a degenerate orbital model, 
g is the Jahn- Teller coupling strength, while in the Hol- 
stein model, it represents the coupling of the lattice to 
the energy of the single orbital present in the model. For 
a strong enough g, an electron becomes localized to the 
lattice site forming a Jahn- Teller (or Holstein) polaron. 
For a given K and g, the magnitude of Q can be min- 
imized to obtain the minimum energy which yields the 
site distortion Qi — grii/K, so that the elastic energy in 
the system becomes EeiasUc = \ Y.i9^^i'^ I ^ ■ 

The Hamiltonian is self-consistently solved within the 
mean-field Hartree approximation : rii-^nn = ni-^{nii) + 
{ni-f)nii — {ni-^){ni^). In the calculations we used a su- 
percell geometry of eleven layers of STO and six layers 
of LAO, with the unit cell doubled along the interface 
in order to study charge disproportionation and polaron 
formation at the interface. The supercell geometry con- 
tained two identical n-type interfaces of the type shown 
in Fig. [T] The ground state energy is given by 

occ 

^ ]^E=-+2Eff'"'V^ (5) 
1 ' 1 ' 

ij ij i 

where the last two terms are the double-counting terms 
and is the number of k points used in the Brillouin 
zone summation. 

Typical parameters for the Hamiltonian are: electron 
hopping -t « -0.15 eV, AT « 5 - 10 eV / A^, g « 2 
eV/A, A = g^/2K w 0.2 eV, and C/ w 2 - 5 eV. We 
have used U = 2 eV in our calculations; the essential ef- 
fect of the Coulomb U in our model is to prevent double 
occupancy of the lattice sites due to the Holstein inter- 
action parameter A. This being the case, a larger U does 
not make any essential change to our result. The di- 
electric constant e can vary over a relatively wide range. 
For example, for bulk SrTiOa, it changes between 25000 
and 300 as temperature changes from 4K to 300Ki^ and 
it decreases sharply with an electric field. The in- 
ternal electric field at the LAO/STO interface, for ex- 
ample, is somewhere between 1-10 mV/A3"— , which 
leads to e ~ 10 — 100 following the non-linear expres- 
sion for e as a function of the electric fieldi^. Our earlier 
density-functional calculations^'^ showed an average di- 
electric constant at the LAO/STO interface to be e ~ 68, 
when lattice relaxation is taken into account, which is 
consistent with the above values. 
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FIG. 2: (Color online) Electronic phases at the polar inter- 
face (schematic). In the 2D phase, electrons are confined 
very close to the interface (one or two monolayers), while 
in the 3D phase, they spread out deeper into the STO bulk. 
Both cases can be metallic (M) or polaronic (P). In the po- 
laronic phase, electrons are localized to individual sites due 
to a strong electron-lattice interaction, which results in an 
insulating phase leading to unpaired electrons on individual 
Ti sites (possible Kondo centers), if the on-site Coulomb U is 
strong enough. 



III. RESULTS AND DISCUSSION 

Fig. [1] shows the typical electron distribution at the 
interface as well as the bare Coulomb potential com- 
ing from the ionic charges (La"'"'^, Al+^, Sr"*"^, Ti+^, 
and O^^), which diverges away from the interface (po- 
lar catastrophe) and the screened potential produced by 
the transfer of the electrons to the interface. Our super- 
cell consisted of two identical n-type Ti02/LaO interface, 
so that exactly half an electron per interface Ti atom is 
available to satisfy the polar catastrophe. As seen from 
the figure, there is no divergence in the screened poten- 
tial, while there is a potential shift due to the accumu- 
lated dipole moment at the interface. 

Competing phases - Our model shows a variety of dif- 
ferent phases depending on the strength of the competing 
interactions, which are illustrated in Fig. [51 If the dielec- 
tric constant e is large, the Coulomb potential of the ions 
(Madelung potential) that confines the electrons to the 
interface becomes weaker, thereby spreading out the elec- 
trons deeper into the SrTiOa side, resulting in a more 3D 
like electron gas. The opposite happens if e is small, so 
that the electrons are confined to just one or two mono- 
layers at the interface, leading to a 2D electron gas. In 
either case, a polaronic phase (2DP or 3DP) can form, if 
the electron-lattice coupling A is strong enough, so that 
electrons are confined to individual sites, giving up the 
kinetic energy they would have gained by delocalization. 
An on-site Coulomb interaction on the Ti site prevents 
the occupation of two such electrons at the same site, so 
that unpaired electrons can form magnetic centers. 
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FIG. 3: (Color online) Phase diagram for the Hubbard- 
Holstein model with half an electron per interface unit 
cell. The star indicates the approximate parameters for the 
LAO/STO interface, which is near several phase boundaries, 
leading to the possibility of phase separation as discussed in 
the text. 



The Hubbard-Holstein model captures the essence 
of the competition between electrostatic and electron- 
lattice coupling and a phase diagram as a function of 
e and A obtained from this model is shown in Fig. |3l 
For the typical parameters for the LAO/STO interface, 
A ~ 0.2 eV and e « 70, the model shows a combined 
3DM -I- 2DM behavior shown by a star mark in the fig- 
ure. In this phase, some electrons are confined to the 
interfacial layers while some are spread deeper into the 
bulk as indicated from the projected densities of states 
(DOS) in Fig. |4l This phase is consistent with the 
density-functional calculations for the interface as well 
as with the experimental studies with very low oxygen 
pressure which shows a mixture of 2D and 3D metallic 
behavior-i^i^i. 

The typical one-electron DOS for the various phases is 
shown in Fig. H) For large e (top left of Fig. H]) , the con- 
fining potential at the interface due to the Madelung term 
is weak, so that the electrons reside both in the interfacial 
and the bulk layers to form a 3D metallic phase. With 
decreasing e, the confining potential becomes stronger, 
leading to an increasing confinement of the electrons to 
the interface layer forming either a 2DM -I- 3DM phase 
or a completely 2DM phase. If the electron-lattice cou- 
pling parameter A is strong enough, the 2D electrons be- 
come localized to individual Ti sites forming a 2D polaron 
phase (bottom right of Fig. S]). 

Charge disproportionation and polaron formation ~ In 
order to examine the instability towards polaron forma- 
tion with the increase of the electron-lattice coupling A, 
we have performed a constrained optimization with re- 
spect to the charge disproportionation using the standard 
Lagrange multiplier technique. The constrained energy 
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FIG. 4; (Color online) Densities-of-states for the interfacial 
layer and the remaining layers (bulk) illustrating the charac- 
teristics of the various phases. Electron states are occupied 
up to the Fermi energy Ef- The letters A and B indicate the 
two inequivalent Ti sites in the interfacial Ti02 plane and 
arrows indicate the two spin states. In the case of the 2DP 
phase for example (bottom right), all electrons are confined 
to the interface layer and occupy the A sites only, one un- 
paired electron per site, while in the 3DM phase (top left), 
the electrons spread out into both the interface sites and sites 
deeper in the bulk. 



is written as 



E{S) = E{{ni}) + T]{nA - ub ~ S), 



(6) 



where the charge disproportionation on the two inequiv- 
alent interface sites, nA — ns, is constrained to the value 
5 and 77 is a Lagrange multiplier. The Lagrange mul- 
tiplier may be written as the partial derivative of the 
constrained energy dE/dS = — ry, so that the constrained 
energy may be obtained by integrating the Lagrange mul- 
tiplier 



r]{S) d6. 



(7) 



The results of the constrained-<5 calculations are shown 
in Figs. [5] and [51 For small A, the minimum energy oc- 
curs a.t S — 0, which indicates that all interfacial Ti sites 
are uniformly occupied. However as A increases, there is 
an instability as the kinetic energy lost by localizing the 
electron at an individual site is overcome by the energy 
gained by the formation of the Holstein polaron. A large 
enough Coulomb repulsion U prevents the accumulation 
of two electrons at the same site. 

Electronic phase separation - The above analysis 
shows charge disproportionation between the two inter- 
facial sites, leading to the A sites which are electron 
rich and B sites which are electron poor. Such electron- 
rich sites may congregate together resulting in electron- 
rich regions and electron-poor regions. In fact, such 
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FIG. 5: Force of constraint 77 as a function of the constrained 
disproportionation parameter 5 (top) and the energy for sev- 
eral values of the electron-lattice coupling strength A (bottom) 
as obtained from the expression Eq. ([7|. Charge dispropor- 
tionation occurs beyond a critical value of A. 
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FIG. 6: Electronic charge in the A and B interfacial sites and 
the rest of the sites (bulk) corresponding to the minima of 
Fig. [5] for various values of A. Beyond a critical value of A, 
a 2DP phase is obtained, where the A sites accommodate an 
electron, leaving the B sites more or less empty. 
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mixed phases are ubiquitous and have been observed 
in many complex oxides including manganite a^^'^^ and 
cobaltates^I and the presence of such inhomogeneous 
phases at the polar interfaces has been suggested to ex- 
plain transport measurements. 

Such inhomogeneous phases can in principle be stud- 
ied by considering very large unit cells, but since the size 
of the individual phase regions could be several nanome- 
ters, the size of such supercells makes such study, either 
from density-functional or from the model Hamiltonian 
studies, computationally prohibitive. To address this is- 
sue, we study the energy of the individual homogeneous 
phases with different electron densities. This however 
makes the individual phases charge non-neutral, making 
their Coulomb energy diverge. But, the phase separated 
system is globally neutral so that the divergences of the 
individual phases cancel out and one need not consider 
these divergences explicitly. The Coulomb interaction 
within and between the phases means that the phases 
would be intermixed in the nanometer scale leading to 
the so-called Coulomb frustrated nanoscale phase sepa- 
ration. What we study here is the tendency of the mate- 
rial to phase separate and not the pattern of the mixed 
phase, which is more complicated requiring a detail bal- 
ance between the Coulomb energy, surface energy of the 
phase boundaries, in addition to the energetics of the 
individual phases. 

To study the electronic phase separation, we varied the 
electron density n of the homogeneous phase and calcu- 
lated the total energy E{n), leaving aside the divergent 
Coulomb energy of the charged phase, which would can- 
cel for the global system as argued above. As a result, 
the expression for the Madelung energy Eq. [5] becomes 

I ' 

T-Lhiad = - ^ Mij{qiqj ~ Sn.iSnj), (8) 

where M is the standard Madelung matrix, qi is the to- 
tal site charge including ionic and electronic of the non- 
neutral system, and Srii = [rii — n^) is the excess charge 
at the i-th site, n,; and being the electronic charge for 
the charged phase and the neutral phase, respectively. 
The above expression includes all Madelung energy ex- 
cept for the Coulomb energy of the excess charge, which 
would diverge with increasing size of the domain. 

In Fig. [7] we have plotted the total energy E{n) with 
respect to the electron concentration n. Note that the 
average concentration for the global system is n = 0.5 
per interface Ti atom. For the stability of a homogeneous 
phase, E{n) must be convex and the figure indicates that 
if the electron-lattice coupling exceeds a critical value 
A > Ac « 0.34, then the system phase separates into an 
electron-rich phase of density ni and an electron-poor 
phase of density n2. A similar conclusion was obtained 
by Kumar et at— from the study of a finite lattice model 
in two dimensions. The chemical potential ^ = dE/dn 
plotted in Fig. [5] shows the phase separated region ob- 
tained from the Maxwell construction. 
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FIG. 7: Energy as a function of electron concentration n of 
the homogeneous phase indicating the existence of phase sep- 
aration if the electron-lattice couphng A exceeds a critical 
value. 
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FIG. 8; Chemical potential as a function of electron density 
at T = and for parameters corresponding to Fig. [7] Phase 
separation occurs beyond the critical electron-lattice coupling 
of Ac. 



To understand the nature of the two phases, we have 
plotted in Fig. [H] the layer distribution of the electronic 
charge for the electron-rich and the electron-deficient 
phases. In the electron-rich phase n = n2 ~ I and the 
electrons occupy the interface Ti sites with occupancy 
nearly equal to one, while in the electron-poor phase 
n = ni, which is nearly zero. These two phases will in- 
termix with the average density of n = 0.5, which is suf- 
ficient to quench the polar catastrophe. Phase 1 may be 
characterized as a low density 3D metallic region, while 
phase 2 is a high-density 2D insulating region. Thus the 
interface region may be partly insulating and partly con- 
ducting leading to a lower number of electrons partici- 
pating in transport as has been seen for example in the 
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FIG. 9: Electron distribution in various layers indicating 
the electron-rich polaronic phase (phase 2) and the electron- 
deficient metallic phase (phase 1) for the case of A = 0.6 eV 
in Fig. El 
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LAO/STO interface. 

Estimation of Tc for phase separation - The tendency 
towards phase separation diminishes with increasing tem- 
perature and disappears beyond a critical value Tc- At a 
finite temperature, one needs to calculate the free energy 
and the chemical potential fJ.{T, n) as a function of tem- 
perature and density. For T = 0, the free energy is / = e 
and /i = d{ne)/dn. For finite T, we compute / and /i 
using a model of non-interacting fermions and assuming 
that only the electronic energy changes with tempera- 
ture. An estimate of T(. may be obtained by using the 
low temperature expansion of the free energy 



(9) 



where a is the coefficient of the Sommerfeld specific heat. 
Furthermore, we approximate the one-particle density of 
states by a semi-circular form 



B K, \ eV being half of the band width, 
potential is then given by 



The chemical 



Ai(r,n)=/i(0) + 



Ef 



6 ^2 



(11) 



where the Fermi energy Ep is a function of the electron 
density n. The results are plotted in Fig. 1101 which shows 
the transition from a mixed phase to a single phase be- 
yond the critical temperature Tc, whose exact magnitude 
depends on the relative strengths of the various interac- 
tions. 

Domain size - So far we have studied the instability of 
the interface to phase separation into two charged phases. 
However, as discussed by many authors, ^^^^i g^^h phase 
separation is frustrated due to the long-range Coulomb 
interaction, leading to a strong nanoscale intermixing be- 
tween the single-phase domains in order to minimize the 
Coulomb energy. The size of the domains in the present 
case may be estimated by minimizing the energy of the 
electron-rich phase. Taking the electron-rich phase to 
have a two-dimensional disk structure, the energy of N 
electrons occupying a single site each in a circle of radius 
R may be evaluated {ttR^ = No?, where a is the lattice 
constant). 

The energy gain due to the Holstein (or JT) polaron 
formation for each electron is —A, while the confinement 
energy of the electron decreases with the domain size 
going as 7 h'^n'^{2mR'^)~^ , where 7 = v'^/tt^ « 0.58, 
ly being the zero of the Bessel function Jo(i^) = 0, and 
the Coulomb energy increases as 2c7^ {3e)~^ , a being 
the surface charge density of the disk. Putting these 
together, the energy per particle E may be expressed 
as 



E = -\ 
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2ma2 



-N- 



2e' 



37r3/2e 



(12) 



D{E) ^ 2n-^B-^VB^ - E^, 



(10) 



the minimization of which yields the optimal domain size 
N = (37/(87r))2/37r3(e^(m*ai)-i)2/3. Taking the relative 
dielectric constant at the STO interface ~ 70, band 
effective mass m* w 10, and lattice constant in units of 
the Bohr radius a/ « 8, the domain size is estimated to 
be TV « 5 or 2R « 10 A. This provides a rough estimate 
of the domain size, which would intermix to minimize 
the long-range Coulomb energy between the electron-rich 
and the electron-poor domains. 



IV. SUMMARY 

To summarize, we studied the Holstein-Hubbard model 
for a polar interface such as LAO/STO in the mean-field 
approximation, which showed the existence of a rich va- 
riety of electronic phases. The phases could be metallic 
or insulating and two or three dimensional depending on 
the strengths of the competing interactions. Our study 
also suggests the formation of phase separation below a 
critical temperature, which makes the transport behavior 
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in these materials quite rich. The phase separation sce- 
nario with mixed metalhc and insulating phase offers a 
possible explanation for the observation of a lower density 
2DEG at the interface than what is necessary to satisfy 
the polar catastrophe. The phase separation is however 
frustrated due to the long-range nature of the Coulomb 



interaction leading to the formation of domains with sizes 
on the scale of several nanometers as has been observed 
in the complex oxides such as the manganites. 

We thank Shi-Jie Chen for stimulating discussions. 
This work was supported by the U. S. Department of 
Energy through Grant No. DE-FG02-00ER45818. 
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